DNA binding analysis of rare variants in homeodomains reveals homeodomain specificity-determining residues.

Homeodomains (HDs) are the second largest class of DNA binding domains (DBDs) among eukaryotic sequence-specific transcription factors (TFs) and are the TF structural class with the largest number of disease-associated mutations in the Human Gene Mutation Database (HGMD). Despite numerous structural studies and large-scale analyses of HD DNA binding specificity, HD-DNA recognition is still not fully understood. Here, we analyze 92 human HD mutants, including disease-associated variants and variants of uncertain significance (VUS), for their effects on DNA binding activity. Many of the variants alter DNA binding affinity and/or specificity. Detailed biochemical analysis and structural modeling identifies 14 previously unknown specificity-determining positions, 5 of which do not contact DNA. The same missense substitution at analogous positions within different HDs often exhibits different effects on DNA binding activity. Variant effect prediction tools perform moderately well in distinguishing variants with altered DNA binding affinity, but poorly in identifying those with altered binding specificity. Our results highlight the need for biochemical assays of TF coding variants and prioritize dozens of variants for further investigations into their pathogenicity and the development of clinical diagnostics and precision therapies.

Temporal dynamics and stoichiometry in human Notch signaling from Notch synaptic complex formation to nuclear entry of the Notch intracellular domain.

Mammalian Notch signaling occurs when the binding of Delta or Jagged to Notch stimulates the proteolytic release of the Notch intracellular domain (NICD), which enters the nucleus to control target gene expression. To determine the temporal dynamics of events associated with Notch signaling under native conditions, we fluorescently tagged Notch and Delta at their endogenous genomic loci and visualized them upon pairing of receiver (Notch) and sender (Delta) cells as a function of time after cell contact. At contact sites, Notch and Delta immediately accumulated at 1:1 stoichiometry in synapses, which resolved by 15-20 min after contact. Synapse formation preceded the entrance of the Notch extracellular domain into the sender cell and accumulation of NICD in the nucleus of the receiver cell, which approached a maximum after ∼45 min and was prevented by chemical and genetic inhibitors of signaling. These findings directly link Notch-Delta synapse dynamics to NICD production with spatiotemporal precision.

Temporal Dynamics and Stoichiometry in Notch Signaling - from Notch Synaptic Complex Formation to NICD Nuclear Entry.

Mammalian Notch signaling occurs when binding of Delta or Jagged to Notch stimulates proteolytic release of the Notch intracellular domain (NICD), which enters the nucleus to regulate target gene expression. To determine the temporal dynamics of events associated with Notch signaling under native conditions, we fluorescently tagged Notch and Delta at their endogenous genomic loci and visualized them upon pairing of receiver (Notch) and sender (Delta) cells as a function of time after cell contact. At contact sites, Notch and Delta immediately accumulated at 1:1 stoichiometry in synapses, which resolved by 15-20 min after contact. Synapse formation preceded entrance of the Notch extracellular domain into the sender cell and accumulation of NICD in the nucleus of the receiver cell, which approached a maximum after ∼45 min and was prevented by chemical and genetic inhibitors of signaling. These findings directly link Notch-Delta synapse dynamics to NICD production with unprecedented spatiotemporal precision.

A spatiotemporal Notch interaction map from plasma membrane to nucleus.

Notch signaling relies on ligand-induced proteolysis of the transmembrane receptor Notch to liberate a nuclear effector that drives cell fate decisions. Upon ligand binding, sequential cleavage of Notch by the transmembrane protease ADAM10 and the intracellular protease γ-secretase releases the Notch intracellular domain (NICD), which translocates to the nucleus and forms a complex that induces target gene transcription. To map the location and timing of the individual steps required for the proteolysis and movement of Notch from the plasma membrane to the nucleus, we used proximity labeling with quantitative, multiplexed mass spectrometry to monitor the interaction partners of endogenous NOTCH2 after ligand stimulation in the presence of a γ-secretase inhibitor and as a function of time after inhibitor removal. Our studies showed that γ-secretase-mediated cleavage of NOTCH2 occurred in an intracellular compartment and that formation of nuclear complexes and recruitment of chromatin-modifying enzymes occurred within 45 min of inhibitor washout. These findings provide a detailed spatiotemporal map tracking the path of Notch from the plasma membrane to the nucleus and identify signaling events that are potential targets for modulating Notch activity.

Lineage-specific control of convergent differentiation by a Forkhead repressor.

During convergent differentiation, multiple developmental lineages produce a highly similar or identical cell type. However, few molecular players that drive convergent differentiation are known. Here, we show that the C. elegans Forkhead transcription factor UNC-130 is required in only one of three convergent lineages that produce the same glial cell type. UNC-130 acts transiently as a repressor in progenitors and newly-born terminal cells to allow the proper specification of cells related by lineage rather than by cell type or function. Specification defects correlate with UNC-130:DNA binding, and UNC-130 can be functionally replaced by its human homolog, the neural crest lineage determinant FoxD3. We propose that, in contrast to terminal selectors that activate cell type-specific transcriptional programs in terminally differentiating cells, UNC-130 acts early and specifically in one convergent lineage to produce a cell type that also arises from molecularly distinct progenitors in other lineages.

Targeted Degradation of the Oncogenic Phosphatase SHP2.

SHP2 is a protein tyrosine phosphatase that plays a critical role in the full activation of the Ras-MAPK pathway upon stimulation of receptor tyrosine kinases, which are frequently amplified or mutationally activated in human cancer. In addition, activating mutations in SHP2 result in developmental disorders and hematologic malignancies. Several allosteric inhibitors have been developed for SHP2 and are currently in clinical trials. Here, we report the development and evaluation of a SHP2 PROTAC created by conjugating RMC-4550 with pomalidomide using a PEG linker. This molecule is highly selective for SHP2, induces degradation of SHP2 in leukemic cells at submicromolar concentrations, inhibits MAPK signaling, and suppresses cancer cell growth. SHP2 PROTACs serve as an alternative strategy for targeting ERK-dependent cancers and are useful tools alongside allosteric inhibitors for dissecting the mechanisms by which SHP2 exerts its oncogenic activity.

IER5, a DNA damage response gene, is required for Notch-mediated induction of squamous cell differentiation.

Notch signaling regulates squamous cell proliferation and differentiation and is frequently disrupted in squamous cell carcinomas, in which Notch is tumor suppressive. Here, we show that conditional activation of Notch in squamous cells activates a context-specific gene expression program through lineage-specific regulatory elements. Among direct Notch target genes are multiple DNA damage response genes, including IER5, which we show is required for Notch-induced differentiation of squamous carcinoma cells and TERT-immortalized keratinocytes. IER5 is epistatic to PPP2R2A, a gene that encodes the PP2A B55α subunit, which we show interacts with IER5 in cells and in purified systems. Thus, Notch and DNA-damage response pathways converge in squamous cells on common genes that promote differentiation, which may serve to eliminate damaged cells from the proliferative pool. We further propose that crosstalk involving Notch and PP2A enables tuning and integration of Notch signaling with other pathways that regulate squamous differentiation.

MAML1-Dependent Notch-Responsive Genes Exhibit Differing Cofactor Requirements for Transcriptional Activation.

Mastermind proteins are required for transcription of Notch target genes, yet the molecular basis for mastermind function remains incompletely understood. Previous work has shown that Notch can induce transcriptional responses by binding to promoters but more often by binding to enhancers, with HES4 and DTX1 as representative mammalian examples of promoter and enhancer responsiveness, respectively. Here, we show that mastermind dependence of the Notch response at these loci is differentially encoded in Jurkat T-cell acute lymphoblastic leukemia (T-ALL) cells. Knockout of Mastermind-like 1 (MAML1) eliminates Notch-responsive activation of both these genes, and reduced target gene expression is accompanied by a decrease in H3K27 acetylation, consistent with the importance of MAML1 for p300 activity. Add-back of MAML1 variants in knockout cells identifies residues 151 to 350 of MAML1 as essential for expression of either Notch-responsive gene. Fusion of the Notch-binding region of MAML1 to the histone acetyltransferase (HAT) domain of p300 rescues expression of HES4 but not DTX1, suggesting that an additional activity of MAML1 is needed for gene induction at a distance. Together, these studies establish the functional importance of the MAML1 region from residues 151 to 350 for Notch-dependent transcriptional induction and reveal differential requirements for MAML1-dependent recruitment activities at different Notch-responsive loci, highlighting the molecular complexity of Notch-stimulated transcription.

Transcriptional Silencers in Drosophila Serve a Dual Role as Transcriptional Enhancers in Alternate Cellular Contexts.

A major challenge in biology is to understand how complex gene expression patterns are encoded in the genome. While transcriptional enhancers have been studied extensively, few transcriptional silencers have been identified, and they remain poorly understood. Here, we used a novel strategy to screen hundreds of sequences for tissue-specific silencer activity in whole Drosophila embryos. Almost all of the transcriptional silencers that we identified were also active enhancers in other cellular contexts. These elements are bound by more transcription factors than non-silencers. A subset of these silencers forms long-range contacts with promoters. Deletion of a silencer caused derepression of its target gene. Our results challenge the common practice of treating enhancers and silencers as separate classes of regulatory elements and suggest the possibility that thousands or more bifunctional CRMs remain to be discovered in Drosophila and 104-105 in humans.

Context and number of noncanonical repeat variable diresidues impede the design of TALE proteins with improved DNA targeting.

Transcription activator-like effector (TALE) proteins have been used extensively for targeted binding of fusion proteins to loci of interest in (epi)genome engineering. Such approaches typically utilize four canonical TALE repeat variable diresidue (RVD) types, corresponding to the identities of two key amino acids, to target each nucleotide. Alternate RVDs with improved specificity are desired. Here, we focused on seven noncanonical RVDs that have been suggested to have improved specificity for their target nucleotides. We used custom protein binding microarrays to characterize the DNA-binding activity of 65 TALEs containing these alternate or corresponding canonical RVDs at multiple positions to ~5,000 unique DNA sequences per protein. We found that none of the noncanonical thymine-targeting RVDs displayed stronger preference for thymine than did the canonical RVD. Of the noncanonical RVDs with putatively improved specificity for guanine, only EN and NH showed greater discrimination of guanine over adenine. This improved specificity, however, comes at a cost: more substitutions of a noncanonical RVD for a canonical RVD generally decreased the protein's DNA-binding activity. Our results highlight the need to investigate RVD-nucleotide specificities in multiple protein contexts and suggest that a balance between canonical and noncanonical RVDs is needed to build TALEs with improved specificity.

Bispecific Forkhead Transcription Factor FoxN3 Recognizes Two Distinct Motifs with Different DNA Shapes.

Transcription factors (TFs) control gene expression by binding DNA recognition sites in genomic regulatory regions. Although most forkhead TFs recognize a canonical forkhead (FKH) motif, RYAAAYA, some forkheads recognize a completely different (FHL) motif, GACGC. Bispecific forkhead proteins recognize both motifs, but the molecular basis for bispecific DNA recognition is not understood. We present co-crystal structures of the FoxN3 DNA binding domain bound to the FKH and FHL sites, respectively. FoxN3 adopts a similar conformation to recognize both motifs, making contacts with different DNA bases using the same amino acids. However, the DNA structure is different in the two complexes. These structures reveal how a single TF binds two unrelated DNA sequences and the importance of DNA shape in the mechanism of bispecific recognition.

Ancient mechanisms for the evolution of the bicoid homeodomain's function in fly development.

The ancient mechanisms that caused developmental gene regulatory networks to diversify among distantly related taxa are not well understood. Here we use ancestral protein reconstruction, biochemical experiments, and developmental assays of transgenic animals carrying reconstructed ancestral genes to investigate how the transcription factor Bicoid (Bcd) evolved its central role in anterior-posterior patterning in flies. We show that most of Bcd's derived functions are attributable to evolutionary changes within its homeodomain (HD) during a phylogenetic interval >140 million years ago. A single substitution from this period (Q50K) accounts almost entirely for the evolution of Bcd's derived DNA specificity in vitro. In transgenic embryos expressing the reconstructed ancestral HD, however, Q50K confers activation of only a few of Bcd's transcriptional targets and yields a very partial rescue of anterior development. Adding a second historical substitution (M54R) confers regulation of additional Bcd targets and further rescues anterior development. These results indicate that two epistatically interacting mutations played a major role in the evolution of Bcd's controlling regulatory role in early development. They also show how ancestral sequence reconstruction can be combined with in vivo characterization of transgenic animals to illuminate the historical mechanisms of developmental evolution.

Direct Promoter Repression by BCL11A Controls the Fetal to Adult Hemoglobin Switch.

Fetal hemoglobin (HbF, α2γ2) level is genetically controlled and modifies severity of adult hemoglobin (HbA, α2ß2) disorders, sickle cell disease, and ß-thalassemia. Common genetic variation affects expression of BCL11A, a regulator of HbF silencing. To uncover how BCL11A supports the developmental switch from γ- to ß- globin, we use a functional assay and protein binding microarray to establish a requirement for a zinc-finger cluster in BCL11A in repression and identify a preferred DNA recognition sequence. This motif appears in embryonic and fetal-expressed globin promoters and is duplicated in γ-globin promoters. The more distal of the duplicated motifs is mutated in individuals with hereditary persistence of HbF. Using the CUT&RUN approach to map protein binding sites in erythroid cells, we demonstrate BCL11A occupancy preferentially at the distal motif, which can be disrupted by editing the promoter. Our findings reveal that direct γ-globin gene promoter repression by BCL11A underlies hemoglobin switching.

Diversification of transcription factor-DNA interactions and the evolution of gene regulatory networks.

Sequence-specific transcription factors (TFs) bind short DNA sequences in the genome to regulate the expression of target genes. In the last decade, numerous technical advances have enabled the determination of the DNA-binding specificities of many of these factors. Large-scale screens of many TFs enabled the creation of databases of TF DNA-binding specificities, typically represented as position weight matrices (PWMs). Although great progress has been made in determining and predicting binding specificities systematically, there are still many surprises to be found when studying a particular TF's interactions with DNA in detail. Paralogous TFs' binding specificities can differ in subtle ways, in a manner that is not immediately apparent from looking at their PWMs. These differences affect gene regulatory outputs and enable TFs to rewire transcriptional networks over evolutionary time. This review discusses recent observations made in the study of TF-DNA interactions that highlight the importance of continued in-depth analysis of TF-DNA interactions and their inherent complexity. This article is categorized under: Biological Mechanisms > Regulatory Biology.

Survey of variation in human transcription factors reveals prevalent DNA binding changes.

Sequencing of exomes and genomes has revealed abundant genetic variation affecting the coding sequences of human transcription factors (TFs), but the consequences of such variation remain largely unexplored. We developed a computational, structure-based approach to evaluate TF variants for their impact on DNA binding activity and used universal protein-binding microarrays to assay sequence-specific DNA binding activity across 41 reference and 117 variant alleles found in individuals of diverse ancestries and families with Mendelian diseases. We found 77 variants in 28 genes that affect DNA binding affinity or specificity and identified thousands of rare alleles likely to alter the DNA binding activity of human sequence-specific TFs. Our results suggest that most individuals have unique repertoires of TF DNA binding activities, which may contribute to phenotypic variation.

Context influences on TALE-DNA binding revealed by quantitative profiling.

Transcription activator-like effector (TALE) proteins recognize DNA using a seemingly simple DNA-binding code, which makes them attractive for use in genome engineering technologies that require precise targeting. Although this code is used successfully to design TALEs to target specific sequences, off-target binding has been observed and is difficult to predict. Here we explore TALE-DNA interactions comprehensively by quantitatively assaying the DNA-binding specificities of 21 representative TALEs to ∼5,000-20,000 unique DNA sequences per protein using custom-designed protein-binding microarrays (PBMs). We find that protein context features exert significant influences on binding. Thus, the canonical recognition code does not fully capture the complexity of TALE-DNA binding. We used the PBM data to develop a computational model, Specificity Inference For TAL-Effector Design (SIFTED), to predict the DNA-binding specificity of any TALE. We provide SIFTED as a publicly available web tool that predicts potential genomic off-target sites for improved TALE design.

Structure of a pseudokinase-domain switch that controls oncogenic activation of Jak kinases.

The V617F mutation in the Jak2 pseudokinase domain causes myeloproliferative neoplasms, and the equivalent mutation in Jak1 (V658F) is found in T-cell leukemias. Crystal structures of wild-type and V658F-mutant human Jak1 pseudokinase reveal a conformational switch that remodels a linker segment encoded by exon 12, which is also a site of mutations in Jak2. This switch is required for V617F-mediated Jak2 activation and possibly for physiologic Jak activation.

DNA-binding specificity changes in the evolution of forkhead transcription factors.

The evolution of transcriptional regulatory networks entails the expansion and diversification of transcription factor (TF) families. The forkhead family of TFs, defined by a highly conserved winged helix DNA-binding domain (DBD), has diverged into dozens of subfamilies in animals, fungi, and related protists. We have used a combination of maximum-likelihood phylogenetic inference and independent, comprehensive functional assays of DNA-binding capacity to explore the evolution of DNA-binding specificity within the forkhead family. We present converging evidence that similar alternative sequence preferences have arisen repeatedly and independently in the course of forkhead evolution. The vast majority of DNA-binding specificity changes we observed are not explained by alterations in the known DNA-contacting amino acid residues conferring specificity for canonical forkhead binding sites. Intriguingly, we have found forkhead DBDs that retain the ability to bind very specifically to two completely distinct DNA sequence motifs. We propose an alternate specificity-determining mechanism whereby conformational rearrangements of the DBD broaden the spectrum of sequence motifs that a TF can recognize. DNA-binding bispecificity suggests a previously undescribed source of modularity and flexibility in gene regulation and may play an important role in the evolution of transcriptional regulatory networks.

Thrombopoietin receptor activation: transmembrane helix dimerization, rotation, and allosteric modulation.

We report how rotational variations in transmembrane (TM) helix interactions participate in the activity states of the thrombopoietin receptor (TpoR), a type 1 cytokine receptor that controls the production of blood platelets. We also explore the mechanism of small-molecule agonists that do not mimic the natural ligand. We show, by a combination of cysteine cross-linking, alanine-scanning mutagenesis, and computational simulations, that the TpoR TM dimerizes strongly and can adopt 3 different stable, rotationally related conformations, which may correspond to specific states of the full-length receptor (active, inactive, and partially active). Thus, our data suggest that signaling and inactive states of the receptor are related by receptor subunit rotations, rather than a simple monomer-dimer transition. Moreover, results from experiments with and without agonists in vitro and in cells allow us to propose a novel allosteric mechanism of action for a class of small molecules, in which they activate TpoR by binding to the TM region and by exploiting the rotational states of the dimeric receptor. Overall, our results support the emerging view of the participation of mutual rotations of the TM domains in cytokine receptor activation.